Interactions and magnetism in graphene boundary states 
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We analyze interaction effects on boundary states of single layer graphene. Near a half filled band, 
both short and long-ranged interactions lead to a fully spin polarized configuration. In addition, the 
band of boundary states acquires a finite dispersion as function of the momentum parallel to the edge, 
induced by the interactions. Away from half filling the wavefunction develops charge correlations 
similar to those in a Wigner crystal, and the spin strongly alternates with the occupation of the 
boundary states. For certain fillings the ground state has a finite linear momentum, leading to the 
formation of persistent currents. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73.43.Lp 
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Introduction. The relativistic low energy properties 
of graphene give rise to many anomalies with respect 
to semiconductor physics and are thus very interesting 
from a fundamental point of view [l], 0]. Furthermore, 
carbon-based nanoelectronics have attracted much inter- 
est as they might complement silicon-based devices 
An important question concerns the magnetic properties 
of graphene. While ideal graphene sheets are far away 
from the ferromagnetic phase transition 0], the occur- 
rence of midgap states, that lead to a peak in the density 
at the Dirac points, can lead to magnetism. 

The most studied example of midgap states are bound- 
ary states 

ESQ- Depending on the kind of edge there 
is a mismatch between the number of A and B sites and 
Lieb theorem (strictly applicable only to Hubbard inter- 
action) guarantees a magnetic ground state at half filling 
with a spin given by S — (N^ — Ng)/2 Q. A prominent 
example are zigzag-edges where the outermost atom cor- 
responds always to the same site. However, boundary 
states are present not only for zigzag edges but for any 
boundary except for pure armchair edges [fj, [9j . The orig- 
inally flat energy band of boundary states is strongly af- 
fected by electron-electron interactions and first principle 



calculations predict magnetic boundaries 

Previous approaches mostly consider the short-ranged 
Hubbard interaction. However, in graphene the screening 
close to half filling is known to be poor [lj] and a com- 
parison with results considering the more realistic case of 
long-ranged Coulomb interactions is therefore desirable. 
For extended systems and finite doping, this is typi cally 
done within the random phase approximation 15|, llq . 
The effect of interaction on the flat band of midgap states 
however is better described by non-perturbative treat- 
ments 



In this work we present a simple effective model which 
allows us to assess the long-ranged Coulomb interaction 
as well as Hubbard interaction. The basic assumption is 
the existence of a filled valence and an empty conduc- 
tion band whose properties do not depend on the fill- 



ing of the boundary states. The model can be justified 
by a) extended states are separated from the degener- 
ate boundary states by an energy gap for any confined 
system, b) minimization of the classical electrostatic en- 
ergy favors the occupancy of the boundary states, c) at 
low energies the density of states is strongly dominated 
by the boundary states. We note that Hartree and ex- 
change interactions between valence band and boundary 
states are taken into account. Finally, the interaction 
between electrons occupying boundary states is treated 
exactly via exact diagonalization of the corresponding 
Hamiltonian. 

Our results are the following: (i) Both interactions give 
similar results close to half filling, however, there are de- 
viations away from half filling, (ii) Our effective model 
conserves all symmetries of the full Hamiltonian and can 
be expressed in terms of low energy states only, (iii) The 
initially flat band of boundary states acquires a finite 
dispersion due to interaction effects, (iv) There is a finite 
regime around half-filling, where the edge is maximally 
polarized, i.e., the edge is magnetic, (v) Strong spin fea- 
tures are predicted for low electron (hole) concentration, 
which should strongly affect the transport properties as 
they may cause spin blockade through graphene islands. 
In the following we first introduce our model and then 
the main results are discussed. Finally we present the 
the conclusions of our work. 

Edge states. While an ideal graphene sheet has a van- 
ishing density of states at half filling, real samples have 
localized states at (or around) zero energy, which are in- 
duced close to edges, impurities, defects or wiggles. As 
a concrete example of a finite set of localized states, we 
choose the boundary states formed at the zigzag edge of 
a carbon nanotube : [B|, 0] 



<t> n = y/l - 4cos 2 (fc„/2)^ [-2cos(fc„/2)] 3 a jkn (1) 

Here, k n = 2wn/N denotes the wave number parallel to 
the edge, N is the number of A-sites at the edge and ajk„ 
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FIG. 1: (Color online) Left: Atomic structure of an un- 
rolled nanotube with zigzag edges, j labels the zigzag lines, 
each containing N different A and B sites, k denotes the 
wavenumber along the edge. Right: Single-particle spec- 
trum for TV = 28, Nb = 9. Red circles denote boundary 
states, black circles extended (bulk) states and blue dashed 
line shows the Dirac cone. Boundary states are labeled by 
integer quantum number m. 



is a plane wave with momentum k n living on the j4-atoms 
of the j-th zigzag line, see Fig. [TJ Boundary states exist 
for all integer numbers n between n min = Int(iV/3 + 1) 
and n max = Int((2iV — l)/3). From now on the boundary 
states are labeled by the integer number m, which we 
call angular momentum, k m — (n m j n + m)2ir/N with 
< m < m max = n ma x ~~ n mm as illustrated in Fig. [TJ 
Nb denotes the number of boundary states. 

We note that the boundary states of Eq. (TIJ are 
the zero energy eigenstates of the original tight bind- 
ing Hamiltonian describing graphene, without any lin- 
earization around the X-points. Thus we also account 
for interaction induced scattering between the two K- 
points, which is particularly important for the bound- 
ary states. Since the system is invariant under exchange 
of the two if-points the energy spectrum is symmetric 
around m sym = m max /2 as illustrated in Fig. [TJ We note 
that boundary states decay fastest for m — m max /2. 

Interactions. We now study the electronic and mag- 
netic properties of the boundary states if either Hubbard 
Vh or Coulomb interaction Vc is included: 

Vh = UJ2^-\) ( n ^-\) W 

Vc = sE^Cni-lHnj-lJ + Vfe (3) 

here i, j denote the lattice sites, m = + the oc- 
cupation of site i with electrons and "f,| the direction 
of the spin. The constant charge density of the posi- 
tive charge background is subtracted from the electronic 
density, which guarantees charge neutrality at half filling 
and establishes electron-hole symmetry around half fill- 
ing. The Coulomb interaction consists of a long-ranged 
part (first part with i ^ j) and the Hubbard part where 



the ratio between both is determined by (e 2 /a)/U (l8j . 

The single particle spectrum of graphene has a sublat- 
tice symmetry. To each states tp v with energy E v < 
there exists a state ip c with energy E c = —E v and vice 
versa. The states are related by ip c (*i) — ^i)ipv( r i) 
with e(rj) = —1 if n is a B-site and e(rj) = 1 otherwise. 
We note that ipv{ r )* ipv ( r ') = ^c(r)*^c( r ')' Therefore we 
can write the completeness relation as 

2^V,(r)' t ^(r')+^V' m (rr^ m (r') = 5 r y , (4) 

v m 

where v denotes the (orbital-)index for the valence band 
and m for the boundary states. 

The interaction between the filled valence band and the 
boundary states can be written as i?vB = ^2 m E m n m , 
where Eq. Q allows to express the potentials E m in 
terms of boundary states only, 

E m — — ^^(nrn\V\nm) + — ~^^(nm\V\mn) . (5) 

n n 

V is either Hubbard Vtf(ri,r 2 ) = U6 Tl , r2 or Coulomb 
interaction Vc"(ri,r 2 ) = e 2 /|ri — r 2 |. Since there is no 
kinetic energy, the effective Hamiltonian for the bound- 
ary states consists only of the potential term Hyb and the 
mutual interaction between electrons occupying bound- 
ary states. We note that the sublattice symmetry al- 
lowed us to account for the interaction with the valence 
band without calculating the finite energy eigenstates. 
In fact, even for a finite energy window around the Dirac 
points one can express all interactions in terms of the 
low energy states within that energy window. The inter- 
action between electrons occupying boundary states is 
treated exactly by numerically diagonalizing the Hamil- 
tonian matrix. Since our effective model conserves the 
symmetries of the full Hamiltonian, we can diagonalizc 
the Hamiltonian matrix in subspaces of given spin, angu- 
lar momentum and particle number denoted as S, S z , M, 
N e . Since the energy is independent of the spin projec- 
tion we set S z = S. Furthermore we note that electron- 
hole symmetry around half filling is also conserved. 

Within our model, the empty band of boundary states 
has zero energy. Since Hy E is diagonal in the basis of 
boundary states the eigenenergy of a single electron oc- 
cupying a boundary state is given by E m ofEq. ©, which 
consists of a Hartree term (first part) and an exchange 
term (second part). One can also define a band structure 
at half filling, where the ground state is maximally spin 
polarized. The energy E m needed to add (or annihilate) 
an electron (or hole) to the ground state at half filling by 
occupying the boundary state \m) is given by 

E m = -^2(nm\V\mn) . (6) 

n 

It only consists of the exchange part. The energy bands 
are depicted in Fig. O 
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FIG. 2: (Color online) Band dispersion E m for empty band 
(see Eq. ©) and half filled band E m (see Eq. ©) for TV = 
28, Nb = 9. Dashed lines show continuous limit (TV — > oo). 



For Hubbard interaction exchange and direct term are 
the same so that E m = —E m with: 



U (l - 4cos 2 (fc m /2)) (1 - 4cos 2 (fc„/2)) 
(1 - 16cos 2 (fc m /2)cos 2 (A:„/2)) 



Here k m = (n m ; n + m)2ir/N as explained below Eq. fTJ). 
These energies rapidly converge for large TV to a band of 
boundary states as shown in Fig. [2j The maximum of 
the band is at m — m sym (k — tt) and the bandwidth 
is given by U(V3/2it - 1/6) w 0.11/7. Close to a K- 
point the dispersion of the boundary states is given by 
zthvphU ~/3t, where k denotes the distance from the K- 
point. 

For Coulomb interaction, exchange and Hartree term 
are different. The exchange interaction is short-ranged 
and converges fast in the limit of a long edge as illus- 
trated in Fig.[2]for the band at half filling. However, the 
Hartree term has a long-ranged contribution that leads 
to a log(TV) divergence of the Hartree contribution to 
the potentials E m . Already the band structure shown 
in Fig. [2] indicates that close to half filling the results 
obtained by either Hubbard or Coulomb interaction will 
resemble each other, while differ results can be expected 
close to an empty or filled band of boundary states. 

Results. As discussed earlier, the effective model for 
the boundary states conserves all symmetries of the full 
Hamiltonian. We note that it is crucial to include the 
Hartree and exchange one-body potentials induced by 
the interaction with the valence band. 

Fig.[3]shows the spin of the ground state (red line) and 
energy gap for spin excitations (black line) as a function 
of the number of electrons occupying boundary states. 
We find that close to half-filling the system is maximally 
polarized, whereas far away from half filling, i.e. close to 
the empty or filled band, the addition of a single electron 
or hole leads to strong alternations of the spin of the 
ground state. The energy gap for spin excitations shows 
an odd-even effect in the particle number and is generally 
bigger close to half filling. 




FIG. 3: (Color online) Spin of the ground state (red) and 
energy gap for spin excitations from the ground state (black) 
as a function of the number of boundary electrons. Upper 
and lower row correspond to TV = 28, TVs = 9 and TV = 
31, TVb = 10 respectively. Energies in units of e 2 /i?(Coulomb) 
and [//TV(Hubbard). 
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FIG. 4: (Color online) TV e -electron ground state energies as 
function of total spin S and angular momentum M (M sym = 
N e m sy m)- The favored spin depends on the angular momen- 
tum. For certain occupation numbers (e.g.TV e = 8) the ground 
state has a finite angular momentum. TV = 28, Nb = 9. 



The underlying physics is determined by the competi- 
tion between the minimization of the interaction energy 
on the one hand and of the single-particle energies on 
the other hand. This competition is most obvious for 
the Hubbard interaction, where only electrons of oppo- 
site spin interact with each other so that the interaction 
between boundary states always favors a spin polarized 
ground state. On the other hand the interaction with 
the filled valence band leads to a dispersion of the band 
of boundary states as illustrated in Fig. [5] and a spin 
unpolarized state allows for a double occupation of ener- 
getically preferred orbitals. 

While close to half filling the spin polarized state is al- 
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ways favored, we find for a nearly empty band of bound- 
ary states that the reduction in potential energy can ex- 
ceed the cost in interaction energy, so that the spin un- 
polarized state is favored. Furthermore, we note that 
there is an odd-even effect in the particle number shown 
in Fig. [31 which is caused by the symmetry of the band 
structure around m sym = m max /2 illustrated in Fig. [2 
which leads to a kind of shell structure in the single- 
particle energies. 

Coulomb and Hubbard interaction strongly resemble 
each other at half filling where the system is charge neu- 
tral and the physics is governed by the short-ranged ex- 
change interaction. Away from half filling, however, the 
structure of the charged edges is dominantly determined 
by charge correlations which are different for Coulomb or 
Hubbard interaction. 

Conclusions and outlook. We have presented a scheme 
to study the electronic structure of midgap states at arbi- 
trary filling, including the effects of interactions. Exam- 
ples are boundary states or states localized near vacan- 
cies, cracks, or wiggles, etc. close to half-filling. In our 
model all interactions are expressed in terms of midgap 
states, and it is exact under the assumption of an inert, 
filled valence band. 

We have applied our model to the localized states at 
zig-zag edges, where the effects of interactions have been 
extensively studied, mostly using short-ranged couplings 
like an onsite Hubbard term and mean field techniques. 
We have not considered in detail effects due to the ex- 
istence of two edges in a graphene ribbon. However, we 
note that the boundary states of different edges are not 
coupled by Hubbard interaction. Treating the kinetic 
coupling perturbatively one then finds that an antiferro- 
magnetic alignment of the two edges is favored, in agree- 
ment with previous work [l(J 0, H3|- We expect the 
same behavior for Coulomb interaction at half filling. 

We found that the ground state close to half filling 
is spin-polarized for both Hubbard and Coulomb inter- 
action, which is in agreement with earlier calculations. 
However, for a low electron (hole) occupation we predict 
strong alternations of the total spin with the number of 
boundary electrons and furthermore we find different be- 
havior for Hubbard or Coulomb interaction. The limits 
can be approximately described as a one dimensional fer- 
romagnet, near half filling, and a Wigner crystal when 
the midgap states are almost empty. 

In both limits, near half filling and an almost empty 
midgap band, interactions lead to an effective one par- 
ticle band which is dispersive, with a bandwidth which 
has a well-defined limit when the length of the edge is 
much larger than the lattice spacing. Hence, the bound- 
ary states acquire a finite velocity, and can contribute to 
the transport properties of the system. We have consid- 
ered clean systems with no disorder. At half filling the 
low energy states have a well-defined valley polarization 
and currents are valley polarized at each edge. All possi- 



ble backscattering processes therefore require intervalley 
scattering, which is usually significantly smaller than in- 
travalley scattering [2l| . 

For certain fillings (see N e — 8 in Fig. [4|), the ground 
state is a doublet with two possible values of the linear 
momentum. If the edge is shaped into a ring, this re- 
sult implies that, for these fillings, the ground state has 
a persistent current of a mesoscopic size, as the current, 
j oc dE/d<p where <fi is an applied flux in a cylindrical 
geometry, scales approximately as Br 1 , where R is the 
length of the edge. As discussed above, the decay of this 
current in samples with disorder is controlled by inter- 
valley scattering. 

Our results also show that the spin of the lowest energy 
state depends strongly on the total linear momentum of 
that state. This situation resembles the dependence of 
the spin on total angular momentum in atoms and quan- 
tum dots with circular or spherical symmetry. 
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